Adiabatic Theory of Nonlinear Evolution of Quantum States 
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We present a general theory for adiabatic evolution of quantum states as governed by the nonlinear 
Schrodinger equation, and provide examples of applications with a nonlinear tunneling model for 
Bose-Einstein condensates. Our theory not only spells out conditions for adiabatic evolution of 
eigenstates, but also characterizes the motion of non-eigenstates which cannot be obtained from 
the former in the absence of the superposition principle. We find that in the adiabatic evolution of 
non-eigenstates, the Aharonov-Anandan phases play the role of classical canonical actions. 
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Adiabatic evolution has been an important method 
of preparation and control of quantum states[[IJ |[ [|. 
The main guidance comes from the adiabatic theorem 
of quantum mechanics B, which dictates that an ini- 
tial nondegenerate eigenstate remains to be an instan- 
taneous eigenstate when the Hamiltonian changes slowly- 
compared to the level spacings. More precisely, the quan- 
tum eigenstate evolves only in its phase, given by the 
time integral of the eigenenergy (known as the dynamical 
phase) and a quantity independent of the time duration 
(known as the geometric phase) . The linearity of quan- 
tum mechanics then immediately allows a precise state- 
ment about the adiabatic evolution of non-eigenstates 
through the superposition principle. 

Our concern here is how the adiabatic theorem gets 
modified in nonlinear evolution of quantum states. Non- 
linearity has been introduced in various forms as possible 
modifications of quantum mechanics on the fundamental 
level ||. Our motivation, however, derives from practi- 
cal applications in current pursuits of adiabatic control 
of Bose Einstein condensates (BECs)0, which can of- 
ten be accurately described by the nonlinear Schrodinger 
equation. Here the nonlinearity stems from a mean field 
treatment of the interactions between atoms. Difficulties 
in theoretical study of adiabatic control of the condensate 
arise not only from the lack of unitarity in the evolution 
of the states but also from the absence of the superposi- 
tion principle. This problem was recently addressed for 
BECs in some specific cases [0, and a similar problem was 
discussed in the past for soliton dynamics 

In this Letter, we present a general adiabatic theory for 
the nonlinear evolution of quantum states (eigenstates or 
noneigenstates). Our study is conducted by transform- 
ing the nonlinear Schrodinger equation into a mathemat- 
ically equivalent classical Hamiltonian problem, where 
nonlinearity is no longer a peculiar issue but rather a 
common character. We can thus apply the adiabatic 
theory for classical systems B |l0[ [ll) to the study of 
adiabatic evolution of quantum states. The eigenstates 
become fixed points in the classical problem, whose adi- 
abatic evolution can be understood from a stability anal- 



ysis of such points. Aharonov-Anandan phases [12| of 
cyclic or quasi-cyclic quantum states play the role of 
canonical action in the classical problem, and are there- 
fore conserved during adiabatic changes of the control 
parameters. In the linear case, i.e., standard quantum 
mechanics, these conserved actions are just the occupa- 
tion probabilities on the eigenenergy levels and the quan- 
tum adiabatic theorem becomes a special case of classical 
adiabatic theorem. Our results are illustrated with a non- 
linear two- level model Jl3|, which describes the tunneling 
of BECs. 

General formalism — The quantum system we con- 
sider follows the nonlinear Schrodinger equation, 

i^^=^(|*),(*|,«)l*(*)). (1) 

where H has been set to one and R denotes all the system 
parameters subject to adiabatic change. For clarity and 
simplicity, we focus in this Letter on quantum systems of 
finite number of levels, that is, = col(V>i, ip2, ••• , ipN) 
and (\ff| = (ipl, -02, •••) ■ We also assume that the 
system is invariant under the global phase transformation 
since it is the case of most interest in physics. The total 
probability (\E r |\f r ) is then conserved and is set to unity 
hereafter. As shown below, it is possible to separate out 
the dynamics of the overall phase, and to describe the 
motion of the remaining variables exactly in terms of an 
(N — 1) dimensional classical Hamiltonian dynamics. 

First, we write |^) = e lX ^ |$), where the overall phase 
A may be chosen to be the phase of any component, say, 
ipff Time evolution of the overall phase is given in terms 
of |$) as derived from Eq.(|l|) 



(2) 



The first term is related to a geometric phase, called the 
Aharonov-Anandan phase []l2f , which will be a key com- 
ponent of our theory. 

Second, without the total phase A, the state |$) be- 
comes a member of the so-called projective Hilbert space 
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and has 2(iV — 1) independent real variables, which can 
be conveniently chosen as 

Q = (l^lM^al 2 ,...,!*^-!! 2 ), (3) 
P = (arg($i),arg($ a ),---,arg($j V _i)). (4) 

They turn out to form a canonical set for an effective 
classical Hamiltonian dynamics as shown below. 

Finally, the dynamics of these variables may be ob- 
tained from the time-dependent variational principle of 
the original nonlinear Schrodinger equation (|lj) with the 
Lagrangian 

C = {V\i-\y)-Hm,{V\,R) 

= ($\i^ t \$)-HmA$\,R)-j t m, (5) 

where Ti. is the total energy of the system and we have 
used the gauge symmetry of the Hamiltonian with re- 
spect to the overall phase. Note that Ti = is 
true only for a linear system; for a nonlinear system, we 
generally have H\Q) = 5H/5(($\). Because the overall 
phase enters in the Lagrangian only through a term of 
total time derivative, the dynamics of the 2(iV — 1) in- 
dependent variables of $) are closed. With the general- 
ized coordinates introduced in Eqs.(3,4) and the relation 
= — Q ■ in virtue of the variational princi- 
ple of the Lagrangian, we find that motions of the gen- 
eralized coordinates are governed by a classical Hamilto- 
nian H d (Q,P,R) = H{\$(JQ,P)),(*(Q,P)\,R), with 
the equations of motion 

dQ = mu dp = _mu 

dt dP ' dt dQ' U 

In this way, we have cast the nonlinear quantum problem 
(|I|) into a formalism of classical dynamics, which will al- 
low us to use some important results on the adiabaticity 
of classical mechanics. 

Eigenstates — We first consider adiabatic evo- 
lution of quantum eigenstates, which can be de- 
fined as extremum states of the system energy 
by H(\&(R)),(&(R)\,R)\&(R)) = Ei(R)\W(R)), 
where the eigenenergies are Lagrange multipliers for im- 
posing the normalization condition of the states. The 
eigenstates correspond to fixed points in the classical dy- 
namics (6) at a given For an elliptic fixed point, 
we expect it to be able to follow adiabatically the con- 
trol parameter provided the latter changes slowly com- 
pared with the fundamental frequencies of periodic orbits 
around the fixed point. The frequencies can be evaluated 
by linearizing Eq. (||) about the fixed point j9[ [ljj, |ll| . In 
standard linear quantum mechanics, these frequencies are 
just the level spacings, so that breakdown of adiabaticity 
occurs by level crossing. In the nonlinear quantum prob- 
lem, the fundamental frequencies are in general different 
from the level spacings, so that adiabaticity can often 



be maintained even if the energy levels cross as demon- 
strated later with a two-level model. 

Nonlinearity in our quantum problem not only makes 
different eigenstates non-orthogonal to each other, but 
also can produce more eigenstates than the dimension N 
of the Hilbert space. Some of these additional eigenstates 
correspond to hyperbolic points in the classical dynamics, 
characterized by complex fundamental frequencies and 
strong sensitivity to small perturbations. We thus ex- 
pect that such quantum eigenstates not to be able to fol- 
low adiabatically the control parameter R. In addition, 
quantum states corresponding to elliptic fixed points can 
collide and annihilate with a hyperbolic point during a 
change of the parameters, leading to the breakdown of 
adiabaticity. 

As an illustration, we consider a nonlinear two-level 
model, 

= (f ° z ~ °* + \ a *) I*) ' W 

where \^) = col(a, b) and a ZjX is Pauli matrix. This 
model was proposed to describe the tunneling of Bose- 
Einstein condensates in optical lattices [l5| or in a 
double- well potentialjjJI. In the model, the parameter 
c characterizes the interaction strength between atoms; 
v is the coupling strength between the two modes. The 
parameter R can be the Bloch wave number |l3|] or energy 
difference between the two wells We are interested 
in the tunneling between the energy levels shown in the 
top panels of Fig|l| when R is increased slowly from the 
far negative end to the far positive end. 

Following the scheme in our general theory, we choose 
the total phase as A = tfb = arg(6) and introduce a pair of 
canonical variables, q = ip a — tfb and p = \a\ 2 . The total 
energy of the system is H = |(*|cr x |*) + f (*|<t z |*) - 
f (*|cr z |*)(*|cr 2 |*). Then we have the equivalent classi- 
cal Hamiltonian as in Eq.(||), 

H cl = v^/p(l-p)cos(q) + |(2p - 1) - - A (2p - If . (8) 

Fig.|l| shows in the top panels the structure of eigenen- 
ergy levels of Eq.(^J) while the phase space orbits of the 
corresponding classical system (||) is shown in the bot- 
tom panels. When c < v, there are only two eigenstates 
and two fixed points. Both fixed points are surrounded 
by periodic orbits and hence elliptic. The frequencies 
of the surrounding orbits approach the fundamental fre- 
quencies of the fixed points as they get closer to the fixed 
points. Since both fx and f% are elliptic, the correspond- 
ing quantum states are expected to be able to follow the 
adiabatically changing R. This has been corroborated by 
our numerical simulations. 

When c > v, the situation becomes very different. In 
the energy band, there are two more eigenstates, forming 
a loop structure; in the phase space portrait, there ap- 
pear two more fixed points with one of them, f^, being 
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FIG. 1: The two top panels show eigen-energies as a function 
of R for two typical cases c < v and c > v. The two bot- 
tom panels show the corresponding phase space portraits at 
a given value of i?(=-0.05). The arrows on the fixed points 
indicate the directions of their movements as R increases. 



hyperbolic. Due to this structure change, the adiabatic 
evolution becomes very different here. First, the eigen- 
state corresponding to /j will not be able to follow the 
adiabatic change of R since is hyperbolic. This has 
been checked by our numerical integration of Eq. (0) . Sec- 
ond, even elliptic point may not able to follow adiabati- 
cally. The fixed point /i can annihilate itself by colliding 
with / 3 as R changes slowly, leading to the breakdown of 
adiabaticity of the tunneling as reported numerically in 
Rcf . . We also notice that there is a level crossing be- 
tween fx and fi at R = 0; however, our calculation shows 
that their fundamental frequencies are ) 2 — l) 1 ^ 2 , not 
equal to zero. This clearly illustrates our statement in the 
general theory that the fundamental frequencies are not 
related to the level spacing in the nonlinear case. 

Cyclic and quasi-cyclic states Compared with 
eigenstates, adiabatic evolution of non-eigenstates is in 
general very complicated as the motions given by Eq. 
@ may be chaotic. We choose to focus on the quan- 
tum states around an elliptic point, where the classical 
orbits are regular. In particular, for the nonlinear two- 
level problem, the classical dynamics is completely in- 
tegrable. The non-eigenstates in such cases are cyclic 
or quasicyclic states, in which the system returns or 
almost returns to its original state after an evolution. 
On the (N — l)-dimensional torus of the regular region, 
we may introduce a set of action-angle variables, with 

i = (h,i 2 ,...,i Nrl ),® = (e 1 ,e 2 ,...e JV _ 1 )@ 0, 0. 

The angular variables change with time with the fre- 
quencies ijj = (u>x, u>2, u>jsr—i) while the actions J" are 
constants. More importantly, according to the classical 
adiabatic theorem [llj , the actions J are adiabatic in- 
variants in the sense that they remain constant even if the 
control parameter R changes (slowly) in time. The ex- 



istence of these adiabatic invariants presents strong con- 
straint on the motion, and guarantees a quantum state 
initially close to an eigenstate (elliptic point) to stay near 
it as the system is changed adiabatically. 

Interestingly, we can attach a physical meaning 
to these adiabatic invariants in the effective classical 
description by making connection to the Aharonov- 
Anandan (AA) phase of the quantum states. The AA 
phase is defined as the time integral of the first term in 
Eq.(|]) for a periodic orbit or a quasi-periodic state (T^|, 



1Aa{R) 



(9) 



We can rewrite it with the canonical variables (Q, P) and 
further with the action-angle variables, 



jaa = f PdQ = i n 



(10) 



where Jl = (o>it, u>2T, ■ ■ ■ ,ujn-xt) and r is time period. 
Therefore, the actions are related to the AA phase jaa, 
which is an observable physical quantity Q . In the spe- 
cial case of N = 2, there is only one independent action, 
so the AA phase is simply Jaa = 27rJ. This simple con- 
nection can be expanded to the general case of TV > 2, 
where one can single out a particular cyclic state that in- 
volves with only one action /„. For this cyclic state, we 
again have the simple relation ^AAn = 2wl n . As a result, 
one can identify these AA phases "fAAn as the adiabatic 
invariants in place of the actions I n . 

How do the above adiabatic invariants connect to the 
familiar notions in the standard linear quantum mechan- 
ics? Consider the time evolution of a general state 
in a linear quantum system for a given R , \^>) = 



Er 



n). where E n, s are the eigenenergies. This 



is a (quasi-) cyclic state with the projective wave function 
given by |*) = Cne'^-^n) + c N \N). Its AA 

phase can be computed with Eq. (||) ; after comparing with 



Eq.fllOD, we immediately find that I n — \c n \ ■ Therefore, 
in the linear quantum mechanics, these adiabatically in- 
variant actions /„ are nothing but the probabilities on the 
energy levels. In this way, we have rederived the adia- 
batic theorem of linear quantum mechanics which states 
that the probability on each energy level is conserved 
in adiabatic processes. Note that this derivation of the 
quantum adiabatic theorem through our effective clas- 
sical description is distinct from the usual semiclassical 
relation discussed in Ref. jl9| . 

The conservation of probabilities on the energy levels 
can be generalized to the case where the system is nonlin- 
ear only in an intermediate range of the parameter R (see 
Fig.^). When the actions (or AA phases) are conserved 
during the entire process, the initial and final probabili- 
ties on the energy levels must remain the same, /„ = I' n , 
because the system is linear at the beginning and end of 
the process so the actions are just the occupation prob- 
abilities. On the other hand, in the intermediate range 
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FIG. 2: Nonlinear tunneling of a system that is nonlinear in 
an intermediate range of the parameter R. Fs and I''s are 
the occupation probabilities on different eigenstates, at the 
beginning and the end of the tunneling process, respectively. 



where the nonlinearity dominates, the probabilities will 
change (even oscillate greatly) since the conserved actions 
are not probabilities on the energy levels. 
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FIG. 3: Left panels: Change of probabilities on the two levels 
with R, which changes with rate a = 0.0001, for three dif- 
ferent cases. The right two panels show how the A A phases 
and the fundamental frequencies change with R in these three 
cases, respectively. 

We now illustrate this important result using our two- 
level model, where, at the two infinite ends with \R\ 3> c, 
the nonlinear term can be ignored and the system is ef- 
fectively linear. For c < v, where all the fixed points 
are elliptic, the fundamental frequency lu for the periodic 
orbit remains finite and the AA phase (action) is con- 
served (see lines (a) in the right panels of Fig|j|). The 
initial and final probabilities on each level are indeed the 



same (FigJ|(a)), although they oscillate in the intermedi- 
ate range of the parameter where the system is nonlinear. 
As the nonlinearity gets strong, the occurrence of tun- 
neling begins to depend on the choice of the initial state. 
In case (b) of Fig.||, where one starts with probability 
/ = 0.1 on level two, tunneling happens; however, in case 
(c) where one starts with probability / = 0.8 on level two, 
there is no tunneling. The difference is whether there is 
collision with the hyperbolic point fs. In case (b), the 
initial noncigenstate falls on a periodic orbit surrounding 
the fixed point f\ , which will later collide with the hyper- 
bolic point f 3 , where the fundamental frequency drops to 
zero and the AA phase has a finite jump (see lines (b) in 
the right panels). The jump height is proportional to the 
tunneling probability. In case (c), the initial state falls 
on a periodic orbit around the fixed point /a, which will 
not collide with f$. 

In conclusion, we have generalized the standard quan- 
tum adiabatic theorem to nonlinear quantum systems by 
transforming them into mathematically equivalent clas- 
sical Hamiltonians. In the classical systems, the eigen- 
states become fixed points and their adiabatic evolutions 
are determined by whether these fixed points are elliptic 
or hyperbolic. Furthermore, we have found that the adi- 
abatic evolutions of non-eigenstates are controlled by AA 
phases, which play the role of classical actions. Within 
the same framework we have also considered the case of 
closed- loop parameter change as pioneered by Berry p0[, 
we find that Hannay's angles|2l| are accumulated at the 
return of the quantum state (details will be published in 
the future). 
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dation, and the International Center of Quantum Struc- 
tures in Beijing. 
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